Signature of Kondo breakdown quantum criticality in optical conductivity 
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We show that the behavior of the finite frequency inter-band transition peak in the optical conduc- 
tivity of the heavy fermions can provide definitive experimental proof of the Kondo breakdown phe- 
nomenon in which the lattice Kondo temperature vanishes at a quantum critical point. Approaching 
such a phase transition from the heavy Fermi liquid side, we find a new cross-over regime where the 
peak position is related to, but is not directly proportional to the lattice Kondo scale. The peak 
position moves to lower energies but remains finite, while the peak value changes non-monotonically 
and it eventually disappears at the quantum critical point. These are unique signatures which 
distinguish a Kondo breakdown transition from a spin density wave driven transition. 

PACS numbers: 71.27.+a, 72.15.Qm, 75.20.Hr, 75.30.Mb 



Introduction. — Developments in material synthesis 
have led to the discovery of numerous rare earth com- 
pounds, the heavy fermions (HFs), which can be tuned 
to a quantum critical point (QCP), separating a magnetic 
ground state from a paramagnetic one [l| , by varying an 
external parameter such as pressure or chemical doping. 
In the quantum critical regime the metallic properties 
are significantly different from what one expects from a 
standard Landau Fermi liquid (LFL), and therefore these 
systems are prototypes to study how strong correlation 
effects give rise to deviations from LFL. 

The early theoretical attempts to describe the QCP 
are based on the possibility that the instability of the 
paramagnetic phase is due to spin density wave (SDW) 
formation 0], and the critical fluctuations are the para- 
magnons. However, in three dimensions these theories 
fail to explain simultaneously the linear temperature (T) 
dependence of the resistivity and the log T dependence of 
the specific heat coefficient observed in experiments 
This has stimulated theorists to construct alternative de- 
scriptions of the quantum criticality [4], among which 
the scenario of Kondo breakdown (KB) has proven to be 
promising [H, [fjj]. In the KB scenario the non-LFL be- 
haviour is due to the presence of a second QCP, in close 
proximity to the magnetic one, where the effective Kondo 
temperature of the lattice (Tk) goes to zero. In this pic- 
ture the critical fluctuations, non-magnetic in origin, are 
associated with the hybridization fluctuation between a 
broad conduction band and a narrow /-electron band. 

In view of the current interest in KB ^ is apt 

to pose the question if a direct experimental signature 
of the KB could distinguish it from other possible routes 
to quantum criticality. The issue is non-trivial partly 
due to the absence of low-T angle-resolved photoemission 
data of the HFs, and partly because standard quantities 
such as resistivity and specific heat do not measure Tk 
directly. For example, the violation of the Wiedemann- 
Franz (WF) law and the divergence of the Gruneisen ra- 
tio (GR) at the QCP have been identified as experimental 



consequences of KB ||. However, a true violation of the 
WF law signifies the absence of electron quasiparticles 
which need not necessarily be due to Tk —> physics. 
An apparent violation of the WF law near a SDW tran- 
sition Q can also take place because of inelastic scat- 
tering. Similarly, the divergence of the GR is expected 
for all kinds of QCPs [l(| ■ Consequently, while KB phe- 
nomenon implies the violation of the WF law and the 
divergence of the GR, the converse does not hold. The 
main purpose of this work is to point out that finite fre- 
quency features in the optical conductivity of the HFs, 
the so-called mid-infrared peak which arise due to inter- 
band optical transitions, can provide a direct and unique 
experimental signature of the KB phenomenon. 

The mid-infrared peak in the optical conductivity of 
the HFs has been well-studied experimentally ll| and 



theoretically [l2[ in the case where the system is far from 
any phase instability. The physical origin of the peak 
is readily understood within a two-band picture of the 
Kondo effect in a lattice, which describes the formation of 
the heavy LFL as the hybridization between a conduction 
c-band and a narrow /-band. The inter-band transition 
takes place between one band with more /-character and 
another with more c-character, and therefore it involves 
an energy-scale at least of the order of the binding energy 
Tk of the composite heavy electrons. Thus the peak 
position is related to Tk, even though the two quantities 
may not always be proportional, as we show below. 

A crucial aspect of the SDW approach is that the para- 
magnetic energy scale Tk remains finite at the QCP. 
Therefore, for such a transition one does not expect the 
mid-infrared peak position and the peak height to change 
significantly while the system is tuned to the QCP by 
varying the external parameter. In contrast, the KB phe- 
nomenon is based upon the possibility of a QCP where 
Tk — > 0. Consequently, a hallmark of KB is a significant 
shift of the inter-band peak position to lower energies, 
as the system is tuned to the QCP. Furthermore, since 
at the KB-QCP the /-band effectively decouples from 
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the c-band, one expects the inter-band feature to gradu- 
ally disappear. These qualitative differences can be use- 
ful to distinguish experimentally a KB transition from 
a SDW instability. This motivates us to study how the 
mid-infrared peak for a heavy LFL varies as the system 
approaches a KB-QCP. 

From the perspective of the charge-dynamics of the /- 
electrons, the KB-QCP has been identified as an orbital 
selective Mott transition, where the /-electrons local- 
ize and decouple from the metallic environment 13 , 3] ■ 
Orbital selective localization is a relevant phenomenon 
in materials with both weakly and strongly correlated 



bands, including classical Mott insulators (like VO2 [15 1 
and V2O3 [IS]), ruthenates Ca2- K Sr :i; Ru04 [13], cobal- 
tates, fullerenes and layered org anic superconductors (for 
a complete list see e.g. Ref. |l8[). KB is likely only one of 
the several possible microscopic mechanisms that leads to 
orbital selective Mott transition. However it is interest- 
ing to conjecture that, besides the HFs, there are other 
cases where the localization can be described in terms of 
a vanishing "effective Kondo temperature". We expect 
that in such systems the study of the optical conductivity 
will be equally germane. 

Model. — The system we study is described by the 
Anderson-Heisenberg model which is given by 



H = — 



(ij),<r i,a 



(1) 



Here t c and tf (t c 3> tf) denote hopping of c- and /- 
electrons respectively, (i, j) are lattice sites (which are, 
say, nearest and next- nearest neighbours) , Eq is the local 
energy level of the /-band, and a is the spin index. Vq 
is the hybridization between the bands, U is the on-site 
Coulomb repulsion in the /-band (U S> tf), and Jr ~ 
(tf) 2 /U is a superexchange spin-spin interaction in the 
/-band. We study the above model in the Kondo limit 
where the /-electrons form local moments. 

The mean field theory is generated by replacing fi a — > 
b\F ia , where (&], &,) describe charged holons and (F- a , 
Fi a ) describe spin-1/2 spinons. The large Coulomb repul- 
sion is taken into account by projecting out the double- 
occupied /-electron states via the constraint F} a Fi a + 
b\bi — 1 (the associated Lagrange multiplier renormalizes 
Eq), At the mean field level this model shows a KB quan- 
tum phase transition at a critical value x c of the parame- 
ter x = Vq/Jh .13]. For x > x c one obtains a heavy LFL 
phase where (hi) ^ and the /-electrons participate in 
band formation, and for x < x c one obtains a phase with 
(hi) = where the /-electrons localize and form a non- 
magnetic uniform spin liquid phase (with (F} a F JrT ) ^ 0). 



As the QCP is approached from the LFL side, the lattice 
Kondo temperature Tjc w n(b) 2 V 2 VQ — > 0, where vq is 
the bare density of states of the c-band. We write the 
spinon bandwidth as at c ~ ((b) 2 tf + Jh), and consider 
a <C 1 as a small parameter of the theory. 

We use linear response theory to study the T = 
optical conductivity of the system as it approaches the 
QCP from the LFL side. We ignore the /-band con- 
tribution to the current operator and write it as J M ~ 

Skcr(^'kek)/jcj C(T Ck ( j, where e^ is the bare dispersion of 
the c-band, and fi denotes spatial direction. This does 
not change our results at a qualitative level, and can 
be formally justified as an expansion in t c /tf. Further- 
more, we write |Vk£k| ~ vp, the Fermi velocity of the c- 
band, and we ignore the vertex corrections to the current- 
current correlator. The optical conductivity ai(tt) is the 
real part of the frequency dependent conductivity tensor, 
which can be written as (for ft > 0, and setting H = 1) 

2e 2 v 2 f° ^ 
o"i(fi) w — —j- / dwVlmGf(w,k)ImGf(w + f7,k), 
J—Q k 

(2) 

where R denotes the retarded function. Within the mean 
field theory the Green's function for the c-electrons is 
G^w.k)" 1 =«-e k + */(2r c )-E«(w,k), £?(w,k) = 
V 2 /[uj - el + i/(2r F )}. V = (b)V is the effective hy- 
bridization between the c-band and the -F-band of the 
spinons, is the dispersion of the spinons whose band- 
width is at c , and r c and Tp mimic the elastic scattering 
due to the presence of impurities. In principle, at c is V 
dependent, but, for simplicity, we treat it as a constant. 

We adopt few simplifications in order to perform the 
calculations analytically: (i) the c-band dispersion is lin- 
earized £k = vp(k — kp), where kp is the Fermi wave- 
vector of the c-band; (ii) the mildly dispersive F-band 



is replaced by a flat level e£ el 



E x [Fig. IHa)] 



Here E x rj avp(kp — kpo), where k F o is the Fermi wave- 
vector of the F-band, and kpo < kp. The latter ap- 
proximation is equivalent to ignoring all non-singular a 
dependence in Eq. ([2]). The microscopic non- universal 
parameter E x , which is the energy necessary for a mo- 
mentum conserving inter-band transition from the Fermi 
surface of the c-band to the F-band, becomes an im- 
portant energy-scale in the limit where the effective hy- 
bridization V is small (see below). We do not expect our 
results to depend on the details of the band-dispersions. 
To demonstrate this explicitly, we also evaluate numeri- 
cally CTi(ri) in Eq. [5]by taking the continuum limit of the 
model. We use parabolic dispersions ek = (k 2 —k 2 F )/ (2m) 
and e£ = (k 2 — k 2 FG )/(2m ), set the energy scale by 
ep = k F /(2m) = 1, and we choose E x = e^/100, 
ro/mo = a = 0.25 We show that at a qualitative level the 
numerical and approximate analytical results match well 
[Figs. [21 E|, and this validates our physical conclusions. 
Results. — Let us for the moment ignore the finite life- 
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FIG. 1: (color online), (a) and e k , are the bare conduc- 
tion and spinon bands (with respective Fermi vectors kp and 
&fo). E x — e1 F is a microscopic energy-scale. E^k are 
the hybridized bands, V is the effective hybridization and 
Q. c = IV is the minimal inter-band transition-energy (corre- 
sponds to momenta for which ek = £k)- (b) For V > E x , the 
optical conductivity peak is at Sip = S7 C . For V < E x the 
fi c -transition is forbidden at T = 0, (since i?ik > is above 
the chemical potential). The peak is at flp — V 2 /E x + E x . 



times of the fermions and write 

ImGf (w, k) = -7r[t4<5(w - E lk ) + ^(w - E 2k )], 




"0 0.05 0.1 0.15 0.2 
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FIG. 2: (color online). Evolution of the optical conduc- 
tivity <Ti(f2) with effective hybridization V ranging from 
0.3 x 10" 2 £ F (left) to 9 x 10~ 2 e F (right), oi(fi) is expressed 
in unit of 00 = e 2 /(mei?). The solid lines (Eq. [2]) are obtained 
using parabolic bands (a = 0.25 and l/r c \m = 10~ 3 ef). The 
dot-dash lines (Eq. [3]) are obtained using a linearized c-band 
and a non-dispersive F-band. In both cases, as V is reduced, 
the peak heights first increase for V 3> E x — ef/100, and 
then diminish for V < E x . For clarity the plots are truncated 
to exclude the Drude feature. 



where u 2 . 
AV 2 ] 1 ' 2 , v } 



[1 + {E x - ek)/D]/2, D = [(e k - E x f + 
= 1— n 2 ,, and Ei, 2 k = [t\t + E x TD]/2 are the 



two hybridized bands. From Eq. ([2]) we get (for il > 0), 



<7i(f2) oc 



n 



duj y u k u k <5(oj 



-n 



k 



£7ik)J(w + O - £ 2k ), 



which is finite provided (a) —51 < £i k < 0, and (b) 
O = i?2k — Eik. Thus, any finite frequency feature in 
the optical conductivity is due to momentum conserv- 
ing inter-band transitions. We note that, for a given 
effective hybridization V, there is a cut-off frequency 
fl c = [E2k — -Eifc] m i n = 2^ below which inter-band tran- 
sitions are not possible. This corresponds to momenta 
where the bare bands e k and e k cross [Fig. Ufa)]. For 
this process u k u k sa 1/4 is maximal, and it reduces for 
inter-band transitions with frequencies larger than Jl c . 

In order to include effects of finite lifetime of the c- 
electrons (we find that finite Tp is of secondary impor- 
tance, and we ignore it in the analytic evaluation), it is 
more convenient to replace the momentum sum in Eq. 
by an energy integral, which can be performed by the 
method of contours for the simplified case of a linearized 
c-band and a non-dispersive F-band. This gives 



«(0) 



2> 



_ n n + i/ Te + - E£(w + n) ' 

(3) 

ne 2 /m, n being the c-electron density. As V 
is reduced we find two different regimes (Figs. [5] and [3]). 



where Oq 



(i) V > E. x . In this regime = Sl c corresponds to 
£i k = E x — V < 0, which implies that the threshold 
frequency for an inter-band transition is Q c [Fig. HJa)]. 
Close to the threshold we get (dot-dashed line in Fig. [2]) 
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fl/2-E x + Z 



n/2- E x - Z 



where Z = {fl\ + i/r 2 ) 1 ' 2 , VL+ = 0(0 2 - Jl 2 ) 1 / 2 /^), 
and 1/t 2 = £1 3 /(4t c S1 2 ). The difference in the posi- 
tion of the peak fl p in txi(il) and the threshold $l c is 
0(l/(Q c r c )), and assuming well defined quasiparticles 
with Q c t c > 1 we get [ Fig. PTJb) and^a)] 



On 



0> 



2V. 



(4) 



Thus, in this regime, which has been discussed in the lit- 



erature 



oc \JTk and the peak position is a di- 
rect measure of the lattice Kondo temperature. For V ^> 
E x , ai(Q p ) oc Slj^Tc/y] 1 / 2 , which implies that the peak 
value increases as V is reduced [Fig. [2] and Fig. EJb)]. 
This behaviour continues until (t c /F) 1 / 2 (^ — E x ) ~ 1, 
after which the peak value decreases with decreasing V. 
We note that, the divergence of <T\ (f2 p ) for infinitely long- 
lived c-electrons is an artefact of a theory where the ef- 
fective hybridization is momentum independent. For the 
case where the momentum dependence of the hybridiza- 
tion is important, we expect <7i(f2 p ) oc 1/V. 
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FIG. 3: (color online), (a) The position f2 p of the inter-band 
peak in oi(fi), and (b) the peak height <ri(fi p ) as a function 
of the effective hybridization V. The triangles in (a) and the 
solid line in (b) are obtained using parabolic bands (Eq.[2]with 
a — 0.25 and 1/t c \f] = 10~ 3 ef). The peak position matches 
the analytical result fi p = ftp (Qp) for V > E x (V < E x ) . 
In (b) the dashed line is obtained using a linearized c-band 
and non-dispersive F-band (Eq. [3)l . 



(ii) V < E x . This new regime, which has not been 
studied earlier, is relevant when the system is close to 
the KB-QCP. We find that ft — il c corresponds to 
_Eik = E x — V > 0, which implies that an inter-band 
transition at this frequency is not possible at T = 0. The 
threshold frequency is il p = V 2 / E x + E x > O c , and this 
corresponds to an inter-band transition from the Fermi 
surface of the lower hybridized band [Fig. [ljb)]. In the 
vicinity of the threshold we get 



( j 1 {n>V 2 /E x + E x ) 



4^V^ 2 ~ n l 



As before, neglecting the finite lifetime, the peak-position 
is at the threshold frequency, and we get [Fig. EJa)] 



n< = V 2 /E x + E x . 



(5) 



Thus, f2 p is no longer a direct measure of Tr-, and in par- 
ticular it stays finite — > E x as Tk — > 0) at the QCP. 
Also, in this regime VL p varies little with V , which may not 
be discernible whenever £l p cannot be very well-resolved. 
Since fi p 7^ fi c , the corresponding matrix element u^v 2 . 
decreases from the maximal value 1/4 as V is reduced. 
As a result the peak- value ai(Q p ) decreases rapidly [Fig. 
[5] and Fig. EJb)], and the peak vanishes at the QCP in- 
dicating the gradual decoupling of the /-band from the 
c-band. Note that if l/r c ^> E x , the inter-band feature 
is eventually masked by the Drude peak. 

Next we discuss the limitations of our results. These 
are based on model calculations and mean field consid- 
erations giving more prominent inter-band feature than 
those observed experimentally. Furthermore, the fluctu- 
ations are non-negligible near the QCP, and the crucial 
question is whether this only smears the feature further 
or washes it out entirely. The issue requires further in- 
vestigation that is beyond the scope of the current paper. 



However, even if the feature is washed out very near the 
QCP, any observation of significant shift of the peak po- 
sition towards lower energy will have demonstrated that 
the physics of these systems is richer than what can be 
captured by SDW theory. 

Conclusion. — We studied the finite-frequency inter- 
band transition peak in the optical conductivity of a 
heavy fermion system near a Kondo breakdown type of 
quantum critical point, where the lattice Kondo temper- 
ature vanishes. As the system approaches the QCP from 
the heavy Fermi liquid side, we find a new cross-over 
regime where the peak position is related to, but is no 
longer a direct measure of the lattice Kondo scale. In 
particular, the peak position moves to lower energies, but 
remains finite at the QCP. On the other hand, the peak 
value changes non-monotonically as Tr- — > 0, and even- 
tually the peak disappears at the QCP indicating the 
decoupling of the /-electrons from the conduction band. 
These features are fingerprints of a Kondo breakdown 
type of QCP, and therefore can be used to distinguish it 
experimentally from a spin density wave type of QCP. 
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